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Abstract 

Clusters of galaxies are the most impressive gravitationally-bound systems in the Universe and 



,^ , its abundance (the cluster mass function) is one important statistics to probe the matter density 

JL. \ parameter (fi^) and the amplitude of density fluctuations (erg). The cluster mass function is 

^ ' usually described in terms of the Press-Schecther (PS) formalism where the primordial density 



fluctuations are assumed to be a Gaussian random field. In previous works we have proposed a 
non-Gaussian analytical extension of the PS approach with basis on the g-power law distribution 
(PL) of the non-extensive kinetic theory. In this paper, by applying the PL distribution to fit 



> 

f~^ ■ the observational mass function data from X-ray highest flux-limited sample (HIFLUGCS) we 



find a strong degeneracy among the cosmic parameters, erg, ^m-, and the q parameter from the 
PL distribution. A joint analysis involving recent observations from baryon acoustic oscillation 
(BAO) peak and Cosmic Microwave Background (CMB) shift parameter is carried out in order 
to break these degeneracy and better constrain the physically relevant parameters. The present 
results suggest that the next generation of cluster surveys will be able to probe the quantities of 
cosmological interest (o"8,Om) and the underlying cluster physics quantified by the g-parameter. 
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I. INTRODUCTION 

The recent astronomical observations are strongly suggesting that the expansion of the 
Universe is speeding up and not slowing thereby revealing the presence of some form of 
repulsive gravity. The basic set of experiments includes: observations from SNe la 1[, cosmic 
microwave background (CMB) temperature anisotropics [2], large scale structure 3|, \^, X- 
ray data from galaxy clustersjsj], age estimates of globular clusters and old high redshift 
galaxies J6|]. In the present cosmic concordance ACDM model the Universe is formed of ~ 
26% matter (baryonic + dark matter) and ~ 74% of a smooth vacuum energy component. 
The thermal CMB component contributes only about 0.01%, however, its angular power 
spectrum of temperature anisotropics encode important information about the structure 
formation process and other cosmic observables. 

On the other hand, the number density of collapsed objects for a given mass at a certain 
time, named the mass or multiplicity function, is a key quantity in the analysis of cosmic 
structures such as clusters of galaxies. The simplest successful approach to analytically de- 
scribe this quantity was developed more than three decades ago by Press and Schechter 7|] 
(hereafter PS). This Gaussian PS formalism is extensively adopted to derive the mass func- 
tion, F(M), of bounded objects in the observed Universeji 

The PS formalism was adopted by Reiprich and Boringer[lll| in their construction of the 
cluster mass function based on the X-ray fiux-limited sample of galaxy clusters (HIFLUGCS) 
selected from the ROSAT All-Sky Survey. As a result, the best fit parameters, Q^n = 0.12 and 
(Tg = 0.96, were obtained in their analysis. These values are, respectively, very low and very 
high when compared with the nowadays independent CMB results [21] thereby leading to some 
skepticism about the usefulness of clusters as sensitive cosmological probes. However, several 
authors have recently claimed that there is no tension between cosmological constraints 



yjipi,y,y,M. 



from CMB and clusters 



12 



131 ]. In particular, Rines and collaborators 13| argued that the 
dynamical determination of cluster masses were overestimated and even modest values of 
velocity segregation between galaxies and dark matter are sufficient to match the mass 
function with the WMAP results. 

In this article we discuss a different possibility. We advocate here that such a discrepancy 
comes out because the matter density fluctuations should be described by an intrinsically 
non-Gaussian random fleld (due to the action of the long-range gravitational interaction) 



and, therefore, the P S appro ach should be somewhat modified. 



,n p^v.p. pape^W y , .n.pi.ed by the T^all. ,-,»„extens.ve stati^cM and kinet.c 
theory 1 171. 1181 . 1191 . |20|, we have proposed a simple extension of the PS analytical formalism. 
Instead the Gaussian function of the original PS approach, it was adopted the power law (PL) 
Tsallis distribution for describing the fiuctuations of the density field. A basic attractive 
feature of the PL distribution is that the resulting model is analytically tractable and the 
standard Press- Schechter formalism is recovered as particular case. In this way, a detailed 
comparison between the two approaches is immediate. 

As we shall see, by applying the PL distribution to fit the observational mass function data 
from X-ray highest fiux-limited sample (HIFLUGCS) we find a strong degeneracy among 
the pair of cosmic parameters {as, fim), and the q parameter from PL distribution. Trough a 



joint analysis involving recent observations from baryon acoustic osci 



and Cosmic Microwave Background (CMB) shift parameter 21 



22 



lation (BAO) peakj^ 



23| the degeneracy is 



broken and the tension between independent determinations is alliviated. 



II. THE PRESS-SCHECHTER APPROACH AND THE Q-POWER LAW 



It is widely known that in the PS original approach the primordial density fiuctuations 
6 = 6p/p for a given mass M is described by a random Gaussian field 
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,j^j. = {5\i) is the mean squared fiuctuation. When the amplitude of the density 
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contrast grows above a critical va^ 

at a given time, can be written as 



ue f^c), a bound object is formed and the fraction F(jv/), 
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while the distribution of bound objects with masses between M and M + dM reads 
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Now, if instead of Gaussian initial fiuctuations, we consider that the amplitudes 
are described by a class of q-parameterized power law distributions, we have the follow 



expression! 
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where the factor B^ is a one-dimensional normahzation constant given by 
a) 5, = (1 - g)^ (V) '^^/^^ , (if < g < 1) 



,(1-9), 



b) 5, = (g-l) 
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(if 1 < g < 2) 



In this framework, the non-extensive multiphcity function of bound objects with masses 
between M and M + dM reads 
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Note that in the hmit q ^ 1 the above PL expressions reduce to the ones of the standard 
Gaussian approach. 



III. POWER LAW DISTRIBUTION AND COSMIC PARAMETERS 

As remarked earher, by performing a x^ statistical procedure with basis on the X-ray 
HIFLUGCS data sample, Reiprich and Boehringer (RB) ll|] determined the statistical con- 
fidence contours for the pair of parameters, erg and ^2^. 

In Figure 1 we show the contours in the Qrn-<^8 plane obtained by using the standard PS 
and PL approaches. The first panel is the RB result. The legend of this panel show the 
parameter q = 1.0, and, as already explained, when the q parameter from the PL distribution 
tends to the unity the Gaussian results are recovered. In this case, the best-fits of the PS 
method are also shown, namely, Qm = 0.12 and as = 0.96. The results from PL distributions 
{q 7^ 1) can be seen in the panels below (g = 1.10, g = 1.15, and g = 1.20, from top to 
bottom). Note that for each panel of Fig{Tl the solid vertical and horizontal lines shows the 
minimum and maximum independent WMAP limits for Vim and erg, respectively 2] • The 
important point here is that the contours using the Gaussian distribution (first panel) does 
not intercept the independent best-fit values from WMAP even at 99% confidence level. In 




FIG. 1: The First panel shows the results from the Gaussian PS approach. The results based on the 
PL distribution appear in the panels below {q = 1.10, q = 1.15, and q = 1.20). The solid vertical 
and horizontal lines show the minimum and maximum WMAP limits for Oj„ and cg, respectively. 
Note that the contours using the Gaussian distribution (the first panel) does not intercept the 
WMAP values while the ones based on the PL distribution are intercepting them for a wide range 
of q values (1 < g < 1.2). 

the other panels, we see that the q free parameter of the PL distribution, however, permit 
the contours to intercept the WMAP resuhs in a wide range of q values. Roughly estimates 
show that the range between q = 1.06 and q = 1.2 is allowed, and, in principle, such a 
degeneracy need to be removed. 

IV. JOINT ANALYSIS AND DISCUSSION 



In principle, the above results suggest that in the non-extensive framework there are 
many possibilities for the theoretical mass function. Still more important, some of them are 
working in the right direction and may help for reconciling the independent estimates of the 
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FIG. 2: Contours on the 0,m-(^8 plane with q = 1.11. The left panel shows poor constraints on the 
cosmological parameters. In the right panel we display the results of the joint analysis with BAO 
and Shift Parameter signatures. 

cosmic parameters. 

In this section we discuss how the parametric 2-dimensional spaces cxs — ^m, ^m — Q 
and (Ts — q can be further constrained by apllying a statistical analyses involving different 
cosmological observations. To this end we consider the current estimates of the baryon 
acoustic oscillations found in the SDSS data |j|, as well as, the shift parameter from WMAP 
observations [2(]. The basic aim is to break the degeneracy between the n^,, erg and q 
parameters in order to better constrain the PL distribution that fits the HIFLUGCS data. 



BAO 



The Baryon Acoustic Oscillations (BAO) in the primordial baryon-photon fluid, leave a 
characteristic signal on the galaxy correlation function, a bump at a scale ~ 100 Mpc, as 
observed by Eisenstein and coworkers ^ . This signature furnishes a standard rule which can 
be used to constrain the following quantity: 
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0.469 ±0.017, 



(6) 



where Ti = H{z)/Hq is the normalized Hubble parameter of the ACDM model, z^, = 0.35 
is a typical redshift of the SDSS sample, and T{z^,) is the dimensionless comoving distance 
to the redshift z^,. This quantity can be used even for more general models which do not 
present a large contribution of dark energy at early times 24 1. 
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FIG. 3: Contours on the ^m-Q plane with ug = 0.81. In the left panel the contour permits almost 
all values for the q parameter. The right panel shows the joint analysis with the BAO and the 
Shift Parameter: the q parameter is restricted in the range 1.0 < q < 1.10, and we have a high 
constraint on the 0,m parameter as well. 

B. CMB shift parameter 



A useful quantity to characterize the position of the CMB power spectrum first peak is 
the shift parameter. For a flat Universe it is given by 

""'■ dz 
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1.70 ±0.03, 



(7) 
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where Zr = 1089 is the recombination redshift and the value for TZ above is calculated from 
the MCMC of the WMAP 3-yr in the standard flat ACDM model. 



C. Results 



Let us know discuss the main results of our statistical analyses. In Fig. [2], by fixing 
q = 1.11, we show the contours on the fi^-crg plane. The left panel shows that we have poor 
constraints on these cosmological parameters. However, when a joint statistical analysis is 
performed using the BAO the shift parameter (right panel) the available space parameter is 
considerably reduced. Note that the contour becomes small and shifted to higher and lower 
values of Qm and ag, respectively. 

By applying the same procedure (now fixing erg = 0.81) and plotting the contours on the 
flm-Q plane, the left panel of Fig. [3] is obtained. One may see that almost all values for the 
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FIG. 4: The left panel shows the contour on the Qm-o's plane when we marginalize over all the 
allowed values of q. In the right panel we apply the BAO and the Shift Parameter signatures. The 
cosmological parameters are in good agreement with the WMAP limits. 

q parameter is allowed, however, a joint analysis with BAO and shift parameter restricts 
severely the q parameter (1 < q < 1.1). Note also that the Q.m parameter becomes tightly 
constrained. 

Finally, we discuss the contours on the Qm-o's plane when we marginalize over all possible 
values of q. The result is shown on the left panel of Fig. HI In the right panel we display 
the result of the joint analysis. It is interesting that by applying BAO and shift parameter 
signature the space parameter is extremely reduced and displaced as occurred in Fig. [2] (see 
also right panel there). We stress that we have marginalized over q and performed a joint 
analysis. The main consequence is that the cosmological parameters are now in agreement 
with the latest WMAP results 



V. CONCLUSIONS 



In the precision era of cosmology it is very important to show the compatibility of the 
cosmological parameters as determined from independent observations. The current stan- 
dard cosmological model, i.e., a flat, accelerating Universe composed of ~ 1/3 of matter 
(baryonic + dark) and ~ 2/3 of a dark energy component in the form of the vacuum en- 
ergy density (A) seems to be fully consistent with a variety of observational data. However, 
some tension has recently been detected between the determinations of the pair {as,fim) 
from galaxy clusters and CMB. In principle, many temporary solutions are possible, but, 



hopefully, only one will survive to future analyses based on the next generation of cluster 
surveys. 

In this article we have discussed an alternative route which seems interesting to be investi- 
gated from a theoretical and observational viewpoint. Initially, inspired by the non-extensive 
kinetic theory, we have extended the original Gaussian distribution of the primordial density 
field|14l. [l5j. By using this PL distribution to fit the observational mass function from the 
HIFLUGCS data we have identified a degeneracy involving the cosmic parameters a^i, and 
Vim with the q parameter from the PL distribution. 

Finally, a joint analysis involving the baryon acoustic oscillation signature and the CMB 
shift parameter have been applied. As a result, we have shown that the non-extensive 
PL distribution may alliviate the tension underlying the independent determinations of the 
cosmic parameters (see right panel of Fig. Hj). 
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